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Abstract 



In this note, we propose a novel application of geoparticle physics, namely using a muon 
tomograph to study ice-filled cleft systems in steep bedrock permafrost. This research could 
significantly improve our understanding of high alpine permafrost in general and climate- 
permafrost induced rockfall in particular. 
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1. Introduction and Summary 

The subject of the proposed research project is the development and application of a muon 
tomograph in order to study ice abundance in high alpine permafrost (figure [T]). To this end 
we want to employ suitable particle detectors to track cosmic-ray muon rates at different 
angles and locations for the purpose of computer-assisted tomographic image reconstruction. 
Differences in count rates from different directions indicate different total masses (density 
lengths) along these directions and can thus be used to create 3D tomographic images using 
computer-assisted tomography (CT). The proposed project constitutes an important new ap- 




Figure 1: Photograph of Matterhorn. Its iconic shape closely resembles our simplified model 
presented below (figure |3| . 

plication of geoparticle physics, i.e., the application of known methods from particle physics 
to study geophysical problems. 
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2. State of research 



2.1. Ice in steep bedrock permafrost 

The stability of many steep mountains is influenced by permafrost (soil or rock material that 
remains at or below zero degrees Celsius for at least two years). Cleft systems that are ice- 
filled and perennially frozen change their stability conditions when warmed [U] or thawed [13] 
and in some cases, this leads to large rock fall events [T^jllf25j . Permafrost thus translates 
climate change into slope stability and affects both the evolution of mountain landscapes 
and natural hazard regimes. 

While considerable progress has been made in understanding the thermal subsurface field in 
mountain topography [15] and in modeling its spatial patterns [16] and transient response 
[2"3] to climate change [21] at diverse scales [2Z] , the distribution of ice-filled cleft systems is 
virtually unknown [13J. 

Currently, two candidate processes are described [13] that can lead to ice-filled cleft systems 
and are investigated in field [TTfB] and laboratory studies [21]. However, their relative im- 
portance and the depth below the surface to which they can act is unknown, despite the fact 
that their presence is the crucial ingredient for climate-permafrost-induced rockfall. Eluci- 
dating the ice content of steep alpine bedrock ridges thus is a key element in understanding 
the mechanisms of ice formation and, especially contributes to anticipating the potential for 
future rock fall from degrading permafrost areas. 

2.2. Geoparticle physics 

Geoparticle physics is the science of probing geological structures by means of cosmic radi- 
ation with methods borrowed from particle physics. There are two subfields defined by the 
particle probe chosen. An important example is neutrino radiography based on cosmic neu- 
trinos. Several neutrino telescopes are under construction, e.g., IceCube under the Antarctic 
ice [3] and the KM3NeT NEMO telescope under the mediteranean sea [18j. Neutrinos also 
feature prominently in some of the observations made at the Pierre Auger Observatory [5|6] . 
Cosmic neutrinos are useful when one wants to probe deep below the Earth's surface or 
explore its core and mantle, due to their extemely large interaction lengths. Therefore they 
are a unique experimental tool to study the Earth's density profile. In contrast to cosmic 
neutrinos, atmospheric muons, i.e., muons produced by the interaction of cosmic rays with 
the atmosphere, can only penetrate through a few kilometers of rock at most, thus they 
can only be used to probe superficial structures such as the internal structure of the Earth's 
crust or mountains. However, to this day, superficial structures are most commonly studied 
by electromagnetic, gravitational or seismologic geophysical techniques which are indirect 
measurements with large intrinsic uncertainties. Muon radiographic or tomographic meth- 
ods present an important alternative tool to study superficial structures such as steep rock 
permafrost. These methods most importantly allow to measure the density of the structure, 
a prominent property that can be readily interpreted in terms of physical composition and 
state. In addition, muon radiography can provide information about the chemical compo- 
sition and thermal state. Especially at sites where the subsurface structure cannot be well 
resolved by conventional techniques due to its inaccessability or heterogeneity, muon detec- 
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tors are suitable alternatives to provide a map of density variations. 

Muon radiography was first proposed by L.W. Alvarez and collaborators in 1970 and em- 
ployed to potentially discover unknown cavities inside the Second Chephren pyramid jl]. A 
similar program is underway in the group of Prof. Schwitters at the University of Texas 
at Austin, where a muon detector is being built for the purpose of exploring the internal 
structure of Mayan pyramids in Belize [12]. Groundbreaking work in muon radiography has 
been done recently by H. Tanaka and collaborators in Japan [28.30,29] . Their observations 
to determine the density distribution of the summit crater of Mt. Asama clearly demonstrate 
the potential of techniques using cosmic-ray muons (figure [2j). 
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Figure 2: (a) Angular distribution of the path length from the topographic map in the top 
region of Mt. Asama. (b) The cosmic-ray muon radiograph is shown in the same angular 
region. Figures taken from [28J. 



2.3. Detector designs 

There are two different types of detectors that can be used to track muon trajectories and 
count rates for the purpose at hand. The mountain environment that the detectors will be 
operated in necessitates special requirements like portability, low or zero electrical power 
consumption and minimal need for maintenance. A full detector setup will consist of many 
detector planes to form a so-called "telescope". 

The detectors serve three primary purposes which are summarized below: 

• Detection of single cosmic-ray muons (events) to enable counting rates at which muons 
hit the target volume. Comparison of rates from different muon directions provides 
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information about absorption rates and thus about the total mass along those muon 
ray paths. 

• Precise measurement of each muon's trajectory to allow a meaningful extrapolation of 
the muon's path through the object to the target volume (which can be assumed to be 
pointlike compared to the object under study). 

• Determination of the muon's energy (velocity) to reject low-energy muons whose large 
scattering would not allow to extrapolate the trajectory to a sufficient degree of accu- 
racy. Alternatively, with certain detector types, it is possible to incorporate a threshold 
which rejects muons whose velocity is below a certain value. 

In the following I describe and compare the advantages and disadvantages of the different 
designs. 

2.3.1. Scintillation detectors (SDs) 

The most common technique for muon detection used in geophysics to date is based on plas- 
tic scintillators strips with wavelength shifting (WLS) fibers frequently used in high energy 
experiments (26] such as MINOS at Fermi lab [1] or OPERA at the Gran Sasso National 
Laboratory [10J. Polysterene scintillator strips are ionization/location detectors that track 
ionization signals which are subsequentally converted into data appropriate for computer 
analysis by employing WLS fibers, photomultiplier tubes (PMT) and front-end electronics. 
The WLS fibers are located in a small groove on the scintillator surface. The scinillation 
light emitted when ionizing particles pass through the scintillator are absorbed by the WLS 
fibers, re-emitted at a slightly different wavelength and transmitted to the PMT tubes at 
both ends with insignificant losses. 

The traditional rectangular design has advantages in being simple, robust, relatively in- 
expensive and well understood by the MINOS group. However, more modern geometries, 
featuring triangular sections, significantly improve the spatial resolution compared to rect- 
angular scintillator designs. With different layering techniques, total sizes of scintillation 
detectors vary from a few m 2 to thousands of m 2 in MINOS and OPERA. 
In order to measure the muon velocity, the UT Muon Maya Group employs a Cherenkov 
detector. This is important to be able to reject low energy muons which feature large scatter- 
ing inside the object under scrutiny, thus deteriorating the theoretical resolution significantly 
(cf. below). However, the size, cost and weight exceed acceptable limits for our type of ap- 
plication. However, there have been significant achievements in reducing size and weight of 
this type of detectors recently. With further improvements, this technology may become a 
feasible option in the future. 

The advantages and disavantages of the scintillator design are listed below. 
Advantages: 

• SDs are commonly used in a wide variety of applications, therefore they are relatively 
inexpensive in production and well understood. 

• Large detector areas can be achieved at low cost. 
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• With moderate detector areas, the power consumption is low enough to be feasible in 
remote environments. The prototype detector of the Vesuvius group in Italy [7] uses 
approx. 70 W of power which can be produced on site, e.g., with solar panels. 

• SDs are extremely robust and insensitive to environmental condtions. Thus they can 
be operated in extreme conditions with minimal on site maintenance. 

• Data acquisition and analysis are very fast. They constitute well studied problems for 
which software packages are readily available. 

Disadvantages: 

• Preliminary estimations of the necessary detector area indicate that, for the purpose 
at hand, the total area should be on the order of magnitude of 100 m 2 to achieve 
reasonable exposure times. While the production of such a detector is possible within 
reasonable bounds of effort and cost, the size makes it difficult to transport and install 
such a system at a remote mountain site. 

• The main disadvantage of SDs compared to the emulsion detectors discussed below is 
their spatial resolution. With the aforementioned setups, angular resolutions of down 
to 10 to 15 mrad in optimal detection conditions seem plausible. This corresponds to 
a spatial resolution of 1 to 1.5 meters at 100 m distance. This is roughly a factor 10 
worse than what can be achieved with nuclear emulsion detectors. 



2.3.2. Nuclear emulsion detectors (EDs) 

Detectors using nuclear emulsion films represent a very old observational technique (origi- 
nally discovered by H. Becquerel, the phenomenon was recognized as being related to the 
ionization of charged particles as early as 1909) in particle physics. Recently, there has been 
renewed interest sparked by technological advances made in Japan during the past few years 
[2811301129] . 

Most importantly, the production technology of the emulsion films (optimization of gel prop- 
erties) and photographic techniques were improved in such a way that the nuclear emulsion 
technique now is feasible for large-scale applications, such as OPERA |22j. Secondly, sig- 
nificant developments in automated emulsion read-out systems nowadays allow for a speedy 
data analysis that proved to be extremely difficult before. The read-out speed is envisaged 
to reach roughly 10 m 2 of film per day per lab in the near future. Such a system consists of 
an automated read-out system, e.g., the Ultra Track Selector (UTS), and off-line software for 
further processing of the huge amount of read-out data. Another important new feature is 
the capability of emulsion- refreshing, i.e., erasing unwanted tracks before use by controlling 
the fading properties of the gel. A detailled account of these improvements can be found in 
[T91I20]. 

Modern nuclear emulsion cloud chambers comprise layers of emulsion film alternating with 
lead plates to achieve a sufficiently high target mass. The emulsion film layers contain mi- 
cro crystals of AgBr, which become developable after a charged particle passes through. 
Therefore a particle trajectory can be tracked after development as a sequence of three di- 
mensionally ordered grains. Emulsion films can contain huge amounts of data stored as 
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photographic images. The real challenge is to convert those pictures into usable data. As 
mentioned above, this challenge has been successfully adressed in recent years. 
Advantages: 

• Highest resolutions of down to 1 mrad can be achieved which correspond to a spatial 
resolution of 10 cm at a 100 m distance. As with SDs, the resolution is naturally 
limited by muon scattering. However, the lead plate layers can effectively be used to 
exclude low-energy muons, thus further improving resolution. 

• EDs are passive devices, i.e., zero power consumption, extremely suitable for the remote 
environment they are supposed to operate in. 

• The size and weight of the EDs are minimal, thus it will be easy to carry them to and 
install them at their intended locations. 

• Due to their small size it is relatively easy to effectively protect the EDs against the 
harsh environmental influences at remote high alpine mountain sites. Moreover they are 
expected to be fully functional in a wide range of temperatures as long as the interlayer 
temperature gradient is not too large, which is not expected to be problematic. 

• Lastly, the costs for production and operation are minimal as well. 

Disadvantages: 

• The read-out and data recovery speed is much slower than conventional systems based 
on PMTs and front-end electronics. However, as mentioned above, there will be sig- 
nificant improvments in the foreseeable future. Moreover, the estimated measurement 
times are on the order of weeks to months, therefore data processing times of days 
to weeks are still acceptable. As part of the OPERA collaboration, LHEP at the 
University of Bern is equipped with the most modern automated emulsion read-out 
microscopes in Europe. Even more capacities will become available shortly. 

In conclusion, the nuclear emulsion detectors seem to be the superior choice of detection 
device for the purpose at hand. In the following we will therefore develop an outline for a 
research program based on this detector technology. 

3. Scientific goals 

The eventual goal of the proposed project is to produce 3D images of high alpine permafrost 
regions inside mountains, using a muon tomograph as described below, in order to determine 
and map their geological composition, with an emphasis on ice content. The importance of 
this was outlined in section (3.1). To this end, I want to formulate the following research 
questions that will serve as a guideline in formulating intermediate goals for the project. 
They will be elaborated further in the next section. 
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Research questions: 

From a geoscience point-of-view: 

• What are the abundance and distribution of ice-filled cleft systems in permafrost of 
steep alpine bedrock ridges? 

• How does climate change influence the slope stability and evolution of mountain land- 
scapes? 

• Can we establish a model of climate-permafrost-induced rock fall in order to identify 
potential future rock fall from degrading permafrost? 

From a physics point-of-view: 

• Can the technology of muon tomography (that is well studied and understood in the 
realm of particle physics) be successfully adapted and applied to a geophysical set- 
ting such as the one proposed herein (with a view toward the difficult environmental 
conditions)? 

• What is the best resolution that can be achieved with nuclear emulsion detectors in 
this framework and what is the maximal thickness of rock that will be accessible under 
realistic conditions? 

• What are the exposure/measurement times necessary to gather reliable data in a high 
alpine setting? 

Intermediate goals: 

• Develop and build a fully functional prototype detector according to specifications and 
durability requirements in order to enable equipment testing under realistic conditions 
and to measure muon fluxes experimentally on site. 

• Identify a suitable mountain ridge for which sufficient topographical data is readily 
available or can be easily obtained to ensure suitable (i.e. small enough) potential 
angular resolutions. 

• Optimize the detector-mountain setup using experimental muon flux and transmission 
data, and Monte Carlo simulations. 

• Development of a computer-assisted tomography software that reconstructs 3D images 
from the data taken by the system of detectors. 



4. Research program 

The main task of the research program is to study ice abundances in steep bedrock per- 
mafrost (with a view towards more general applications of studying the internal structure 
of mountains). To achieve this, we propose to employ a suitable emulsion detector system 
to track cosmic-ray muon rates at different angles (and different locations) for the purpose 
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of 3D tomographic image reconstruction. This is possible because differences in count rates 
from different angles indicate different density lengths along those directions which allows to 
detect internal structures (such as ice-filled cleft systems with a typical density of almost f g 
cm~ 3 ) embedded in the surrounding rock (with an assumed standard density of 2.5 g cm -3 ) 
and enables computer-assisted tomographic reconstruction. A simplified model sketch of the 
setup is drawn in figure [3] 




Figure 3: Sketch of the mountain-detector setup. The depicted idealized internal structure 
is assumed to be spherical (despite the fact that cleft systems will be planar features). 

The research program has to be structured and devided into several main stages of 
research and development. 

At first it will be necessary to identify a suitable object of study, i.e., a potentially interesting 
mountain ridge or steep rock face that fulfills our requirements concerning size and shape. Its 
size should not exceed a few to several hundred meters of rock. Moreover, we need sufficient 
topographical data at an accuracy of < 5 cm to be able to achieve precision measurements 
in the vicinity of 1 mrad angular resolution. The technique presented herein is especially 
suitable for thin rock layers, i.e., little rock thickness. The reason is that, for muon energies 
above 100 GeV, the muon flux has a power law dependence on the energy, with exponent 
-2.7, and thus falls off rather steeply. Therefore the observation of muons traversing thick 
layers of rock becomes increasingly difficult. The energy loss of relativistic particles in matter 
can be estimated to be roughly 2.5 MeV g _1 cm 2 . Assuming a standard rock density of 2.5 
g cm -3 , the energy loss is a little more than 0.5 TeV km~ l . The integrated muon flux as a 
function of rock thickness is plotted in figure [4] (taken from |30j.) 

Furthermore, it will be necessary to adapt the nuclear emulsion detectors to our purposes, 
most importantly to install them in a weatherproof casing and to make sure they can with- 
stand the extreme temperature differences and environmental conditions on site. The next 
step is to determine the optimal specifications and locations for the detectors. This includes 
a simulation of the experimental apparatus and estimation of muon fluxes. To this end, one 
has to do feasibility studies starting with a Monte Carlo simulation to estimate the various 
muon fluxes and transmission rates. The optimization process comprises gradual changes of 
detector design and locations and will lead to the final optimal design and locations of the 
emulsion detectors. For the Monte Carlo simulations, a wide variety of software packages 
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Figure 4: Integrated flux of cosmic ray muons at various zenith angles pentrating through a 
given thickness of rock with a density of 2.5g/cm 3 . 

(e.g., GEANT4 |2j) is readily available and has been tested in many (sometimes similar) 
applications. A digital elevation map of the mountain will be the geometrical input to the 
Monte Carlo simulation. In a first step, under some hypothetical assumption of the internal 
structure, one generates a large number of tracks which run backwards from the detector 
and intersect the digital elevation map and the hypothetical substructure to be 'revealed'. 
For each tracjectory, the different layers of matter and their density are recorded. This data 
subsequently serves as input for GEANT4 with the aim of determining muon transmission 
fluxes. This should be supplemented by actual experimental muon measurements on site, 
helping to understand the experimental problems and to validate the Monte Carlo simula- 
tions. At the end of this stage one should be able to produce first radiographical images that 
could be compared to conventional geophysical techniques. 

Thirdly, simplified model calculations of the mountain-detector system are necessary to es- 
timate the detector size (effective area) and exposure time. A rough estimate based on the 
simplified model depicted in figure [3] suggests that exposure times on the order of a few 
weeks are realistic. This will likely be improved once the exact detector configuration has 
been determined and allows for realistic simulations and testing. 

Ultimately, the goal is to install several detectors at optimized locations to be able to recon- 
struct 3D tomographical images of the internal structures. This final step will take a strong 
effort of software development but, fortunately, tomographic methods are widely used in 
medical and particle physics and are well studied and documented. 

Acknowledgements 

I am indebted to Stefan Gruber for suggesting this project, numerous discussions and proof- 
reading an earlier version of this note. I would also like to thank Antonio Ereditato, Sonia 
Paban and Roy Schwitters for discussions and comments. 



9 



References 

[1] P. Adamson, G. Crone, L. Jenner, R. Nichol, R. Saakyan, C. Smith, J. Thomas, M. Ko- 
rdosky, K. Lang, P. Vahle, A. Belias, T. Nicholls, G. Pearce, D. Petyt, M. Barker, 
A. Cabrera, J. Hartnell, P. S. Miyagawa, N. Tagg, A. Weber, E. Falk Harris, P. G. 
Harris, R. Morse, P. Symes, D. Michael, P. J. Litchfield, R. Lee, and S. Boyd. The 
MINOS calibration detector. Nuclear Instruments and Methods in Physics Research A, 
556:119-133, January 2006. 

[2] S. Agostinelli et al. GEANT4: A simulation toolkit. Nucl. Instrum. Meth., A506:250- 
303, 2003. 

[3] J. et. al. Ahrens. Sensitivity of the IceCube detector to astrophysical sources of high 
energy muon neutrinos. Astroparticle Physics, 20:507-532, February 2004. 

[4] L. W. Alvarez, J. A. Anderson, F. El Bedwei, J. Burkhard, A. Fakhry, A. Girgis, 
A. Goneid, F. Hassan, D. Iverson, G. Lynch, Z. Miligy, A. Hilmy Moussa, Mohammed- 
Sharkawi, and L. Yazolino. Search for Hidden Chambers in the Pyramids. Science, 
167:832-839, February 1970. 

[5] L. Anchordoqui, T. Han, D. Hooper, and S. Sarkar. Exotic neutrino interactions at the 
Pierre Auger Observatory. Astroparticle Physics, 25:14-32, February 2006. 

[6] Luis A. Anchordoqui, Haim Goldberg, Dan Hooper, Subir Sarkar, and Andrew Tay- 
lor. Predictions for the cosmogenic neutrino flux in light of new data from the pierre 
auger observatory. Physical Review D (Particles, Fields, Gravitation, and Cosmology), 
76(12):123008, 2007. 

[7] F. et. al. Beauducel. Muon radiography of volcanoes and the challenge at Mt. Vesuvius. 

[8] J. Beutel, S. Gruber, A. Hasler, R. Lim, A. Meier, C. Plessl, I. Talzi, L. Thiele, 
C. Tschudin, M. Woehrle, et al. PermaDAQ: A scientific instrument for precision sens- 
ing and data recovery in environmental extremes. In Proc. 7th Intl Conf. Information 
Processing Sensor Networks (IPSN09), page to appear. ACM Press, New York, 2009. 

[9] M.C.R. Davies, O. Hamza, and C. Harris. The effect of rise in mean annual tempera- 
ture on the stability of rock slopes containing ice-filled discontinuities. Permafrost and 
Periglacial Processes, 12(1), 2001. 

[10] M. Dracos. The OPERA experiment. Physics of Atomic Nuclei, 67:1092-1096, June 
2004. 

[11] M. Geertsema, J.J. Clague, J.W. Schwab, and S.G. Evans. An overview of recent large 
catastrophic landslides in northern British Columbia, Canada. Engineering geology, 
83(1-3):120-143, 2006. 

[12] UT Maya Muon Group. A detector for muon tomography - technical report. 



10 



[13] S. Gruber and W. Haeberli. Permafrost in steep bedrock slopes and its temperature- 
related destabilization following climate change. Journal of Geophysical Research, 
112(F2), 2007. 

[14] S. Gruber, M. Hoelzle, and W. Haeberli. Permafrost thaw and destabilization of Alpine 
rock walls in the hot summer of 2003. Geophysical Research Letters, 31(13):L13504, 
2004. 

[15] S. Gruber, M. Hoelzle, and W. Haeberli. Rock-wall temperatures in the Alps: mod- 
elling their topographic distribution and regional differences. Permafrost and Periglacial 
Processes, 15(3):299-307, 2004. 

[16] S. Gruber, L. King, T. Kohl, T. Herz, W. Haeberli, and M. Hoelzle. Interpretation 
of geothermal profiles perturbed by topography: The Alpine permafrost boreholes at 
Stockhorn Plateau, Switzerland. Permafrost and Periglacial Processes, 15(4):349-357, 
2004. 

[17] A. Hasler, I. Talzi, J. Beutel, C. Tschudin, and S. Gruber. Wireless Sensor Networks 
in Permafrost Research: Concept, Requirements, Implementation, and Challenges. In 
Proc. 9th Intl Conf. on Permafrost (NICOP 2008), volume 1, pages 669-674, 2008. 

[18] Ulrich F. Katz. KM3NeT: Towards a km**3 Mediterranean neutrino telescope. Nucl. 
Instrum. Meth., A567:45 7-461, 2006. 

[19] K. Kodama et al. Observation of tau-neutrino interactions. Phys. Lett., B504:218-224, 
2001. 

[20] K. Kodama et al. Detection and analysis of tau neutrino interactions in DONUT 
emulsion target. Nucl. Instrum. Meth., A493:45-66, 2002. 

[21] J.B. Murton, R. Peterson, and J.C. Ozouf. Bedrock fracture by ice segregation in cold 
regions. Science, 314(5802):1127-1129, 2006. 

[22] T. Nakamura, A. Ariga, T. Ban, T. Fukuda, T. Fukuda, T. Fujioka, T. Furukawa, 
K. Hamada, H. Hayashi, S. Hiramatsu, K. Hoshino, J. Kawada, N. Koike, M. Komatsu, 
H. Matsuoka, S. Miyamoto, K. Miyanishi, M. Miyanishi, K. Morishima, H. Nada, N. Na- 
ganawa, T. Nakano, K. Narita, M. Natsume, K. Niwa, N. Nonaka, B. D. Park, O. Sato, 
S. Takahashi, T. Toshito, T. Uetake, M. Nakamura, K. Kuwabara, S. Nishiyama, 
Y. Nonoyama, and K. Kodama. The OPERA film: New nuclear emulsion for large-scale, 
high-precision experiments. Nuclear Instruments and Methods in Physics Research A, 
556:80-86, January 2006. 

[23] J. Noetzli and S. Gruber. Transient thermal effects in Alpine permafrost. The 
Cryosphere, 3(l):85-99, 2009. 

[24] J. Noetzli, S. Gruber, T. Kohl, N. Salzmann, and W. Haeberli. Three-dimensional 
distribution and evolution of permafrost temperatures in idealized high-mountain to- 
pography. Journal of Geophysical Research, 112(F2), 2007. 



11 



[25] J. Noetzli, M. Hoelzle, and W. Haeberli. Mountain permafrost and recent Alpine rock- 
fall events: a GIS-based approach to determine critical factors. In 8th International 
Conference on Permafrost, Proceedings, edited by: Phillips, M., Springman, S., and 
Arenson, L, volume 2, pages 827-832. 

[26] A. Pla-Dalmau, A. D. Bross, and K. L. Mellott. Low-cost extruded plastic scintillator. 

Nuclear Instruments and Methods in Physics Research A, 466:482-491, July 2001. 

[27] D. Riseborough, N. Shiklomanov, B. Etzelmuller, S. Gruber, and S. Marchenko. Recent 
advances in permafrost modelling. Permafrost and Periglacial Processes, 19(2):137-156, 
2008. 

[28] H. K. M. Tanaka, T. Nakano, S. Takahashi, J. Yoshida, and K. Niwa. Development 
of an emulsion imaging system for cosmic-ray muon radiography to explore the inter- 
nal structure of a volcano, Mt. Asama. Nuclear Instruments and Methods in Physics 
Research A, 575:489-497, June 2007. 

[29] H. K. M. Tanaka, T. Nakano, S. Takahashi, J. Yoshida, H. Ohshima, T. Maekawa, 
H. Watanabe, and K. Niwa. Imaging the conduit size of the dome with cosmic-ray 
muons: The structure beneath Showa-Shinzan Lava Dome, Japan. Geophysical Research 
Letters, 34:22311-+, November 2007. 

[30] H. K. M. Tanaka, T. Nakano, S. Takahashi, J. Yoshida, M. Takeo, J. Oikawa, T. Ohmi- 
nato, Y. Aoki, E. Koyama, H. Tsuji, and K. Niwa. High resolution imaging in the 
inhomogeneous crust with cosmic-ray muon radiography: The density structure below 
the volcanic crater floor of Mt. Asama, Japan. Earth and Planetary Science Letters, 
263:104-113, November 2007. 



12 



